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Abstract 

The adsorption of highly charged flexible polycations and polyanions on a charged cylindrical 
substrate is investigated by means of Monte Carlo (MC) simulations. A detailed structural study, 
including monomer and fluid charge distributions, is provided. The influence of a short range 
attraction between the polycations and the negatively charged substrate is also considered. We 
demonstrate that the building up of multilayer structures is highly prohibited mainly due to the 
high entropy penalty stemming from the low dimensionality of the substrate at strong curvature. 
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I. INTRODUCTION 



The adsorption of charged polymer chains [polyelectrolytes (PEs)] on oppositely charged 
surfaces and particles has been the subject of intense theoretical and experimental research 
work. Due to the strong mutual electrostatic attraction, the resulting complex is highly 
stable and can exhibit new and interesting physico-chemical properties. For instance, such a 
process is employed to stabilize charged colloidal suspensions. Another more elaborated PE 
adsorption technique, consists of successive deposition of polycations (PAs) and polyanions 
(PAs) leading to well controlled PE multilayer structures. Such materials present enormous 
technological applications and this deposition method has been especially exploited with 
planar 1 and spherical substrates. In principle this technique could also be applied to (ultra- 
thin) cylindrical substrates. However due to the very low dimensionality there, a quantitative 
experimental characterization of the PE multilayering is more difficult. 

On the theoretical side, there exist a few analytical works about PE multilayers on charged 
planar surfaces,—^ but no theory is available for thin cylindrical substrates. Recent MC 
simulations were devoted for the PE multilayering on planar— and spherical charged surfaces. 
The aim of this paper is to fill this gap and to investigate by means of MC simulations the 
structure of PCs and PAs adsorbed onto a charged cylindrical surface. The relevant role of 
the entropy is especially pointed out. 

Our paper is organized as follows: Sec. HU is devoted to the description of our MC 
simulation model. The measured quantities are specified in Sec. IIHI The PC adsorption is 
studied in Sec. |TV| Then the adsorption of both PC and PA chains is addressed in Sec. |V| 
Finally, Sec. I VII provides brief concluding remarks. 

II. SIMULATION MODEL 

The setup of the system under consideration is similar to those recently investigated with 
spherical 7 and planar^ substrates. Within the framework of the primitive model we consider 
a PE solution near a charged hard rod with an implicit solvent (water) of relative dielectric 
permittivity e r ~ 80. This fixed charged rod (whose axis is located at r = \J x 2 + y 2 = 0) 
is characterized by a negative linear bare charge density — Aoe, where e is the (positive) 
elementary charge and Aq > is the number of charges per unit length. Electroneutrality is 
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always ensured by the presence of explicit monovalent (Z c = 1) rod's counterions of diameter 
a. PE chains (N + PCs and AL PAs) are made up of N m monovalent monomers (Z m = 1) 
of diameter a. Hence, all micro ions are monovalent: Z = Z c = Z m = 1 and have the same 
diameter a. For the sake of simplicity, we only consider here symmetrical complexes where 
PC and PA chains have the same length and carry the same charge in absolute value. Here 
the radius r rod of the charged rod is also a/2. To keep further simplicity in our model, 
we also assume the same linear charge density for the rod — Aoe = —le/a as the flexible 
polyelectrolytes. Hence the only differences between the rod and the PE chains are the 
chain stiffness and and the chain length, where both of them are infinite for the rod. 

All these particles making up the system are confined in a R x L cylindrical box. Periodic 
conditions are applied in the z direction, whereas hard coaxial cylinders are present at 
r = a/2 (location of the cylindrical charged wall) and at r — R (location of an uncharged 
cylindrical wall). To avoid the appearance of image charges,^ we assume that on both parts 
of the charged rod (at r = a/2) the dielectric constants are the same. 

The total energy of interaction of the system can be written as 



where the first (single) sum stems from the interaction between an ion i (located at r = r^) 
and the charged plate, and the second (double) sum stems from the pair interaction between 
ions % and j with R^ = |Rj — Rj| where Rj = Xie x + yie y + ZiB z . All these contributions to 
Utot in Eq. (JTJ) are described in detail below. 

Excluded volume interactions are modeled via a hardcore potential defined as follows 




(1) 



[Uhs(Rij) + U cou i(Rij) + U FE NE{Rij) + U LJ (Rij)] , 



i,i<j 





for the microion-microion one, except for the monomer-monomer one,-^ and 



oo, for n < r rod + a/2 
Uhs° d \zi) = < oo, for n > R- a/2 

0, for r rod + a/2 < n < R - a/2 



(3) 
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for the rod-microion one. Hence the minimal radius, ro, of closest approach between the 
rod-axis and a microion is ro = r ro d + a/2 = a. 

The electrostatic energy of interaction between two ions % and j reads 

U C oul [Rij j 



knT 



(4) 

iB-L JTLij 

where +(-) applies to charges of the same (opposite) sign, and Is = e 2 /(47ie e r k B T) is the 
Bjerrum length corresponding to the distance at which two monovalent ions interact with 
k B T. The electrostatic energy of interaction between an ion % and the (uniformly) charged 
rod reads 

±2l B X \n(ri/r rod ) (5) 



U {rod) (r- 

U coul l'« 



where +(-) applies to positively (negatively) charged ions. An appropriate and efficient mod- 
ified Lekner sum was utilized to compute the electrostatic interactions with periodicity in one 
direction.— To link our simulation parameters to experimental units and room temperature 
(T = 298K) we choose a = 4.25 A leading to the Bjerrum length of water lg = 1.68a = 7.14 
A and to the so-called Manning condensation paramete r 12 ! 13 £m = Ib/cl = 1.68. 

The polyelectrolyte chain connectivity is modeled by employing a standard FENE poten- 
tial in good solvent (see, e.g., Ref.— ), which reads 



U 



FENE 



(Rij) 



--nRn In 



R J 



for < Rn 



(6) 



oo, 



for > Rq 

with k = 27ksT ' j 'a 2 , Rq = 1.5a, and % and j stand for two consecutive monomers along a 
chain. The excluded volume interaction between chain monomers is taken into account via 
a purely repulsive Lennard- Jones (LJ) potential given by 

12 



U LJ (R, 



4e 



0. 



ij 



'j 



+ e, for R^ < 2 1 / 6 



for R^ > 2 1 / 6 a 



(7) 



where e = ksT. These parameter values lead to an equilibrium bond length I = 0.98a. 

An important interaction in PE multilayering is the non-electrostatic short ranged at- 
traction, Uf^fa , between the cylindrical macroion and the PC chain. To include this kind 



TABLE I: Model simulation parameters with some fixed values. 



Parameters 



T = 298K room temperature 

— Aoe = — 40e/L = — e/a macroion linear charge density 

Z = 1 microion valence 

a = 4.25 A microion diameter 

I = 0.98a bond length 

r rod = a/2 rod radius 

r o = r rod + a/2 = a radius of closest approach 

Ib = 1.68a = 7.14 A Bjerrum length 



R = 40a (\/x 2 + y 2 )-box radius 

L = 40a z-box length 

N + number of PCs 

iV_ number of PAs 

N PE = N + + total number of PEs 

N m = 20 number of monomers per chain 

Xvdw strength of the specific VDW attraction 



of interaction, we choose without loss of generality a (microscopic site-site) van der Waals 
(VDW) potential of interaction between the cylindrical macroion and a PC monomer that 
is given by 

Uvd°w ( r *) = -tXvdw ( _l — ) for Ti > r , (8) 
\ ' i ' o \ a j 

where Xvdw is a positive dimensionless parameter describing the strength of this attraction. 
Thereby, at contact (i.e., r = r ), the magnitude of the attraction is Xvdw^ — XvdwksT which 
is, in fact, the relevant characteristic of this potential. 1 - Since it is not straightforward to 
directly link this strength of adsorption to experimental values, we choose Xvdw = 5 as 
considered previously.—^ 

All the simulation parameters are gathered in Table The set of simulated systems can 
be found in Table |H] The equilibrium properties of our model system were obtained by us- 
ing standard canonical MC simulations following the Metropolis scheme.— 1 ^ Single-particle 
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TABLE II: System parameters. The number of counterions (cations and anions) ensuring the 
overall electroneutrality of the system is not indicated. The chain length N m = 20 is always the 
same as well as the parameters of the cylindrical macroion. 



System N PE N + A_ 

A 6 6 

B 12 6 6 

C 24 12 12 

D 48 24 24 



moves were considered with an acceptance ratio of 30% for the monomers and 50% for the 
counterions. Typically, about 10 5 to 10 6 MC steps per particle were required for equilibra- 
tion, and no less than 5 x 10 5 subsequent MC steps were used to perform measurements. 
To improve the computational efficiency, we omitted the presence of PE counterions when 
N + = iV_ so that the system is still globally electroneutral. We have systematically checked 
for N + = AL = 6 (system B) that the (average) PE configurations (especially the monomer 
distribution) are qualitatively the same of those where PE counterions are explicitly taken 
into account, as it should be. 

III. MEASURED QUANTITIES 

We briefly describe the different observables that are going to be measured. In order to 
characterize the PE adsorption, we compute the monomer density n±(r) that is normalized 
as follows 

R-a/2 

n±{r-)2nrLdr = N±N m (9) 

where +(— ) applies to PCs (PAs). This quantity is of special interest to characterize the 
degree of ordering in the vicinity of the cylindrical macroion surface. 

The total number of accumulated monomers N±(r) within a distance r from the axis 
(located at r = 0) of the cylindrical macroion is given by 

N ± (r) = / n±(r')2irr'Ldr' (10) 

Jr 
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where +(— ) applies to PCs (PAs). This observable will be addressed in the study of PC 
adsorption. (Sec. IIV|) . 

Another relevant quantity is the global net fluid charge A(r) which is defined as follows 



where n + stand for the density of all the positive (negative) microions (i.e., monomers 
and counterions) . Thus, A(r) corresponds to the (reduced) net fluid charge per unit length 
(omitting the bare macroion linear-charge — Ao) within a distance r from the axis of the 
charged cylindrical wall. At the uncharged wall, electroneutrality imposes X(r = R — a/2) = 
Ao- By simple application of the Gauss' law, [A(r) — Ao] fr is directly proportional to the 
mean electric field at r. Therefore A(r) can measure the screening strength of the cylindrical 
macroion-charge by the neighboring solute charged species. 

IV. ADSORPTION OF POLYCATIONS 

In this part, we study the adsorption of PC chains (system A) for the two couplings Xvdw = 
and Xvdw = 5. The case Xvdw = 0, corresponding to a purely electrostatic regime, is well 
understood for "normal spherical" electrolyte ions near a charged rod . 1211 ^ 1 ? Nevertheless, 
the situation is much more delicate for PE adsorption, and that has not been the subject 
of numerical simulations.— This is a decisive step to elucidate the even more complex PE 
multilayer structures where additionally PAs are also present. 

Here, where iV_ = (i.e., no polyanions), global electroneutrality is ensured by the 
presence of explicit PC's counterions (i.e., monovalent anions) and the macroion-rod's coun- 
terions (i.e., monovalent cations). 

The profiles of the monomer density n + (r) are depicted in Fig. ^ Near contact (i.e., 
r — r ~ 0), the density n + {r) at Xvdw = 5 is about one order of magnitude larger than at 
Xvdw = 0. To further estimate the degree of adsorption of PC monomers we have plotted the 
fraction N + (z)/(N + N m ) of adsorbed monomers in Fig. El At a radial distance r — ro = a 
(corresponding to a width of two monomers), almost all the the monomers (more than 99%) 
are adsorbed for Xvdw — 5 against only ~ 68% for Xvdw = 0. This strong observed PC 
adsorption at Xvdw = is of course due to the high valence of the polycation (here N m = 20) 
and the intense attractive electric field near the thin rod. On the other hand, at high 
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FIG. 1: Profiles of PC monomer-density ?i_i_(r) at different \ydw 

couplings (system A). 
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FIG. 2: Fraction N + (z)/(N + N m ) of adsorbed PC monomers at different Xvdw couplings (system 
A). 

enough Xvdw and (even) without rod-monomer electrostatic interaction [i.e., for a neutral 
rod (A = 0)], the PC adsorption would be very similar to that obtained at finite A . 

The (global) net fluid charge A(r) is reported in Fig. EJ In all cases we detect a macroion- 
surface charge reversal (i.e., A(r)/A > 1). The position r = r* at which A(r*) gets 
its maximal value decreases with Xvdw, due to the Xw2w- en hanced adsorption of the PCs. 
Concomitantly, this overcharging increases with Xvdw, since the (extra) gain in energy by 
macroion-monomer VDW interactions can better overcome (the higher Xvdw) the cost of 
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FIG. 3: Net fluid charge A(r) at different x-vdw couplings (system A). 



the self-energy stemming from the adsorbed excess charged More quantitatively, we have 
A(r* — r ~ l.la)/Ao ~ 2.96 at Xvdw = 5 against A(r* — r ~ 3.4a)/Ao ~ 2.67 at Xvdw — 0. 
Note that the maximal value of charge reversal of (6 x 20-40) /40 = 200% (i.e., A(r*) /A = 3) 
allowed by the total charge of PCs can not be exactly reached due to a slight accumulation 
of microanions. However at Xvdw — 5, the presence of microanions near the rod is very 
marginal, leading in practice to an overcharging that is fully governed by the PCs. Again 
this giant overcharging found at any Xvdw is due to the high (intrinsic) electric field generated 
by the macroion rod. Indeed, for the same chain parameter, it was found in planar geometry 
(with a surface charge density — a e = — 0.165C/m 2 ) that only a charge reversal of about 
70% [i. e., a(r*)/cr = 1.7] could be reached at Xvdw = 5.— For cylindrical substrates, it turns 
out that the equivalent surface charge density [here — A e x L/(27rr rod L) « — 0.282C/m 2 ] 
can be much higher in experimental conditions. Those (locally) overcharged states should 
be the driving force for the building of subsequent PE bilayers when PA chains are added. 

Typical equilibrium configurations can be seen in Fig. 0] The qualitative difference 
between Xvdw = [Fig. Ufa)] and Xvdw = 5 [Fig. Efb)] is remarkable. Without additional 
VDW attraction {Xvdw — 0) the adsorption is more diffuse than at Xvdw — 5, where in the 
latter situation the r-ffuctuation is very weak. Interestingly, one can observe that the PC 
chains tend to adopt a helical structure. This kind of conformation is the non-trivial result of 
two antagonistic forces: (i) electrostatic correlations and (ii) entropy. To clarify this point, 
one has to consider two typical length scales: 
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FIG. 4: Typical equilibrium configurations for PC chains adsorbed onto an oppositely charged 
cylindrical macroion (system A). The little counterions are omitted for clarity. The outer green 
cylinder is a guide for the eye. (a) Xvdw = (b) Xvdw = 5. 

• If the cylindrical macroion is larger or about the PE chain (i.e., L > N m l as it is 
presently the case), then the ground stated corresponds to straight PE chains sticked 
to the surface of the rod. At finite temperature one obtains PC chains spiraling around 
as the result of the coupled effects of entropy and electrostatic correlations. 

• If L <C N m l, then the ground state (already) corresponds to perfectly ordered heli- 
cal structures (so as to overscreen the rod) which locally resemble a one-dimensional 
Wigner crystal.— At finite temperature and strong Coulomb coupling (as it is the case 
with highly charged PEs and macroion-substrates) , this high helical ordering persists 
at shorter range. 

Those reasonings above should at least hold when the equivalent surface charge densities of 
the PE and that of the rod are similar, as it is the case in this study. 
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FIG. 5: Profiles of monomer density n±{r) for oppositely charged polyelectrolytes (systems B — D). 
(a) Xvdw = 0. (b) Xvdw = 5. 

V. ADSORPTION OF POLYCATIONS AND POLYANIONS 

We now consider the systems B — D where additionally PA chains are present, so that 
we have a (globally) neutral polyelectrolyte complex (i.e., N + = AL = Npe/2). Global 
electroneutrality is ensured by the counterions of the cylindrical macro ion as usual. We 
stress the fact that this process is fully reversible for the parameters under investigation. In 
particular, we checked that the same final equilibrium configuration is obtained either by (i) 
starting from system with isolated PCs and then adding PAs or (ii) starting directly with 
the oppositely charged PEs. 

A. Monomer distribution 

The profiles of the monomer density n±(r) at Xvdw = and Xvdw = 5 are depicted in Fig. 
Efa) and Fig. Efb), respectively. The corresponding microstructures are sketched in Fig. H3 

1. Case of zero-Xvdw 

At Xvdw = 0, a comparison between systems A (see Fig. and B indicates that the 
adsorption of PC monomers is weaker when additional PAs are present. This effect was 
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FIG. 6: Typical equilibrium configurations for the adsorption of oppositely charged PE chains 

(systems B — D) onto a cylindrical macroion. The polycations are in white and the polyanions in 

red. The little ions are omitted for clarity. The outer green cylinder is a guide for the eye. (a) 

Xvdw = 0, N PE = 12 (system B) (b) Xvdw = 5, N PE = 12 (system B) (c) Xvdw = 0, N PE = 24 

(system C) (d) Xvdw = 5, N PE = 24 (system C) (e) Xvdw = 0, N PE = 48 (system D) (f) Xvdw = 5, 

12 

N PE = 48 (system D). 



already observed with spherical and planar— substrates, and the same mechanisms apply 
here to cylindrical substrates. More explicitly, the PC tends to build up a globular state 
(reminiscent of the classical bulk PE collapse) by getting complexed to the PA chain, as 
perfectly illustrated in Fig. HJa). Thereby, the mean number of nearest monomer- neighbors 
gets higher which is both (i) entropically and (ii) energetically (at least from the PE complex 
viewpoint) favorable. It turns out that for macroion-rods with large curvature, this PC 
desorption is counter-intuitively much higher than in the planar case,— although a much 
weaker charge reversal (without PAs) was found for the latter. Indeed, we just saw that 
the overcharging by PCs alone can be quite huge (of the order of 200 % even at Xvdw = 
- see Fig. EJ), and thus, one would naively expect a new wave of charge reversal in system 
B (where N + = N_ = 6) of the same order (see also Sec. IV Bl for a more quantitative 
description) which in turn should lead to a (very) stable bilayering. But Figs. Efa) andEfa) 
demonstrate that it is not the case. Furthermore, an inspection of the density profile n_(r) 
of the PA monomers reveals that there is no peak for Npe = 12. Those interesting findings 
allow us to already draw a partial conclusion: 

• At large enough rod-curvature and small enough PE concentration, even if a giant 
overcharging would occur in the absence of PAs, it is not sufficient to promote bilay- 
ering. 

Hence, this spectacular PA-induced PC desorption is due to entropic effects. More pre- 
cisely, the adsorption onto quasi one-dimensional objects is extremely costly from an entropy 
viewpoint, and if the number of available dipoles (i.e., PC-PA ion pairs) is too small then 
ion-dipole correlations are not strong enough to overcome the "entropic barrier" . 

Upon increasing Np% (i.e., the PE concentration) one sees that the density n + (r) at con- 
tact increases [also identifiable on Figs. EJa), (c) and (e)] and gets gradually less sensitive 
to Npe, showing that a "steady" regime should be reached at large enough Npe- Concomi- 
tantly, by increasing Np&, Fig. EK a ) shows that the adsorption of PA monomers increases 
[also identifiable on Figs. Eta), (c) and (e)]. This iVp E -enhanced adsorption of monomers is 
essentially governed by two concomitant mechanisms: 

• At higher Npe, the global ion-dipole correlation gets enhanced (due to an increase of 
PC-PA ion pair concentration) which favors a stronger monomer adsorption. 
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• Concomitantly, at higher Npe, where the available volume becomes lower, the entropy 
loss undergone upon adsorption gets reduced, which also favors a stronger monomer 
adsorption. 

Nevertheless, only a very broadened peak in n_(r) (at Npe = 48) is seen and no second peak 
in n + (r) appears, indicating a weakly stable bilayer. This feature strongly contrasts with 
our previous findings for planar substrates^ where stable trilayers were reported at Xvdw = 
(also with N m = 20). Therefore, we can state: 

• PE multilayering in the electrostatic regime is highly inhibited at large rod-curvature, 
mainly due to entropic effects. It is to say that the (conformational) PE entropy along 
the rod is much smaller than that in the bulk and also smaller than that on a planar 
surface. 

2. Case of finite \vdw 

In general all our mechanisms previously described at Xvdw = are also present here 
at Xvdw — 5, but one has now to take into account the possible dominance of the VDW 
interaction near the surface of the cylindrical macroion. 

At Xvdw = 5, Fig. E(b) shows that the first PC layer is always very stable (despite of 
the presence of PAs) due to the the sufficiently strong rod-PC VDW attraction. In parallel, 
n + (r) and in general the ordering increase with Npe, as it should be. At large Npe = 48, 
one can detect a second broadened peak in n + {r) which could be the signature of the onset 
of a third PE layer. However, the height of this peak is still significantly smaller than 
the corresponding decaying PA density n_(r) in this region. A visual inspection of Fig. 
Efb), (d) and (f) confirms all those features. Hence, even with the help of a specific VDW 
attractive interaction, it seems that the formation of stable structures beyond bilayers is 
highly inhibited at strong confinement. This again proves how important is the role of 
entropy in PE multilayering with low dimensional substrates. 

B. Charge distribution 

The net fluid charge A(r) is reported in Fig. [7| In all cases, the cylindrical macroion gets 
(locally) overcharged. Surprisingly, beyond the PC overcharging layer the fluid charge does 
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FIG. 7: Net fluid charge A(r) at different Xvdw couplings (systems B — D). 

not reenter below unity (at least for r < 6a), proving that there are no charge oscillations 
in the vicinity of the rigid rod. Even with an overcharging of A(r*)/Ao = 3.5 for Xvdw = 5 
and Npe = 48, no charge oscillation is generated. This is, nevertheless, consistent with our 
previous results of Fig. 03 where no multilayering was observed. This phenomenon is due to 
the fact that our PE concentration is still (very) small even at Npe = 48, which is also the 
experimental situation in the process of PE multilayering where a rinsing is applied at each 
successive layer deposition^. It is well known that at high finite electrolyte concentration, 
charge oscillations take place near a charged cylindrical wall . 1 ^ 2 ? 124 

Our results clearly show that the degree of overcharging is systematically larger at high 
Xvdw — 5, as also observed without PAs (see Fig. |HJ). The values of r* are also systematically 
shifted to smaller ones as one increases Xvdw However, the position of the peak (r* — r ~ a) 
for Xvdw = 5 remains nearly independent of N PE in contrast to what happened at Xvdw = 0, 
where r* is systematically shifted to smaller radial distances upon increasing Npe- All those 
observations are consistent with our mechanisms presented so far. 

Nonetheless, as soon as oppositely charged polyions can interact, there is a subtle interplay 
between clustering and the lateral correlations of polyions that governs the degree of over- 
charging near the macroion surface, as it was already demonstrated in Refs.^. Additionally, 
at strong rod curvature, entropic effects become also considerable. Consequently, the physics 
of overcharging involved in the adsorption of oppositely charged PEs onto one-dimensional 
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substrates is highly complicated. 



VI. CONCLUDING REMARKS 

In summary, we have carried out MC simulations to elucidate the adsorption of flexible 
highly charged polycations and polyanions onto a charged cylindrical substrate. Within 
the dilute regime, we considered the influence of PE concentration. In order to enhance 
the possible formation of PE multilayers we have also considered an extra non-electrostatic 
short-range attraction (characterized here by Xvdw) between the cylindrical macroion surface 
and the polycations. 

As far as the PC adsorption is concerned (in the absence of PAs), we demonstrated that 
huge macroion charge reversal occurs even in a purely electrostatic regime with Xvdw = 0. 
By adding exactly the same amount of PAs, we surprisingly observe a (relatively) marginal 
overcharging which is due to (i) PC-PA clustering and (ii) above all to entropic effects. 

At higher number of PEs, our results show that true bilayering (i.e.; flat and dense PE 
layers) can only occur at finite Xvdw, m contrast to what was recently found with planar 
substrates.— Even at finite Xvdw, our simulation data demonstrate that stable multilayering 
(beyond bilayering) is hard to reach at large macroion rod-curvature, due to the high entropy 
loss there. This latter in turn inhibits the appearance of charge oscillations. 

However, there must be a certain regime of curvature, typically when the rod radius is 
larger or about the chain size (i.e., r rod > N m l), where such entropy effects become less 
relevant and would then lead to similar behaviors to those occurring with planar substrates. 
In this respect, a future work could include the effect of rod-curvature on the formation 
of PE multilayers. We hope that this contribution will trigger new theoretical as well as 
experimental investigations. 
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